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ABSTRACT 

We investigate the physical and chemical evolution of population II stars with initial 
masses in the range 6.5-8 M Q , which undergo an off centre carbon ignition under 
partially degenerate conditions, followed by a series of thermal pulses, and supported 
energetically by a CNO burning shell, above a O-Ne degenerate core. In agreement 
with the results by other research groups, we find that the O-Ne core is formed via the 
formation of a convective flame that proceeds to the centre of the star. The evolution 
which follows is strongly determined by the description of the mass loss mechanism. 
Use of the traditional formalism with the super-wind phase favours a long evolution 
with many thermal pulses, and the achievement of an advanced nucleosynthesis, due 
the large temperatures reached by the bottom of the external mantle. Use of a mass 
loss recipe with a strong dependence on the luminosity favours an early consumption 
of the stellar envelope, so that the extent of the nucleosynthesis, and thus the chemical 
composition of the ejecta, is less extreme. The implications for the multiple populations 
in globular clusters are discussed. If the "extreme" populations present in the most 
massive clusters are a result of direct formation from the super-AGB ejecta, their 
abundances may constitute a powerful way of calibrating the mass loss rate of this 
phase. This calibration will also provide informations on the fraction of super-AGBs 
exploding as single e-capture supernova, leaving a neutron star remnant in the cluster. 

Key words: Stars: abundances - Stars: AGB and post-AGB 



1 INTRODUCTION 

The deep spectroscopic surveys of Globular Clusters (GC) 
stars performed in the last decades have revealed star- 
to-star differences, that trace well defined abundance pat- 
terns, involving all the "light" elements, up to Aluminium 
l|Norris et al.l fl98ll ) . The discovery that stars belonging to 
the main sequence or the sub-giant branch show the same 
pattern s of abundance s present in the red giants of the same 
cluster l|Carrettall2006r ) . indicated that some self-enrichment 
mechanism has been active in GCs, and that a new, sec- 
ond, stellar generation (SG) formed, from the ashes of the 
evolution of older objects, belonging to a first generation 
(FG). A key-point towards the understanding of the pro- 
cesses behind these observational results is the large fraction 
of stars with the anomalous chemistry, that hardly drops 
below 50% of the total num ber of stars examined , and ex- 
ceeds 80% in some clusters (|Carretta et al J l2009al lbh . The 
Oxygen-Sodium anticorrelation is by far the most studied, 
and has been confirmed to exist in practically all the GCs in- 
vestigated, although the slope of the O-Na relationship, and 
the lowest oxygen abundances detected, vary considerably 
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from cluster to cluster. In many clusters a Mg-Al anticor- 
relation is also detected, althoug h the extent of the magne - 
sium depletion is still debated l|Cohen fc Melendej 12005). 
Unlike the O-Na trend, the Mg-Al relationship is not of 
straight interpretation, because not only the slope (if any), 
but also the maximum magnesium and the minimum alu- 
minium mass fractions, which are commonly interpreted as 
the abundances of the primordial popul ation, change from 
cluster to cluster (jCarretta et al.lr2009bl h Photometric sig- 
natures are generally less striking than the spectroscopic 
evidence, apart from the presence of very anomalous "blue" 
main sequences (bMS) in tw o among the most massive clus- 
ters ( ui Cen and NGC 2808 jBedin et al1l2004l; IPiotto et ail 
l2005l : iD'Antona et~al1 120051 : IPiotto et al1l2007D and by the 
presence of very extended horizontal branches and very lu- 
minous RR Lyrae in two metal rich clusters -NGC 6388 
and NGC 6441. Both these features indicate the presence of 
a very heliu m rich population (Y ^ 0.36-0.38) in t hese four 
clusters (e.g. lNorrij[2004 ICaloi fc D Antonall2007l ). 

Three distinct scenarios have been proposed so far for 
the progenitors of the SG, each of which must face some 
difficulties in reconciling the theoretical predictions with 
the observational ev idence: a) winds from massive AGBs 
l|Ventura et al.l l200ll). b) ejecta from fastly rotating mas- 
sive stars ( Decressin et al.l 120071 ) , 



and c) massive binaries 
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l|de Mink et al.|[2009l ). All these polluters together, plus the 
possible contribution of s t ellar collisions, have been pro- 
posed bv lSills fe Glebbeekl l|2010h . In this work wc focus on 
the first hypothesis, alrea d y out lined on qualitative grounds 
by I Ventura fc D'Antonal (2OO8I ) and IVentura fc D'Antonal 
(2009), and fixed on a most robust framework by two semi- 
nal papers, computing first the possible hydrodynamical and 
N-body evolution, and further exploring the chem ical evo- 
lution of protoclusters l|D'Ercole et al.ll200Sl . [2010) ). In these 
latter works, it is suggested that the SG forms in a cooling 
flow established in the GC core by the gas ejected at low 
velocity by massive AGBs. This process is restricted to a 
narrow time interval 10 8 yr, thus allowing only the con- 
tribution of stars with masses 5M©, after which star 
formations stops, due to the onset of SNIa explosions. The 
loss of stars from the outskirts of the cluster helps diminish- 
ing the FG/SG stars number ratio, because this mechanism 
hardl y involves, in a first p hase, the co re born SG objects . 

ID'Ercole etHI (|2008l ), following iPumo et alj l|200Sl) . 
also suggested that in the most massive clusters a premi- 
nent role in the star forming process is played by super- AGB 
(SAGB) stars, defined as objects with initial masses in the 
range (9-11M0 j]] that undergo off-center carbon ignition in 
partially degenerate conditions, and end-up their evolution 
as O-Ne white dwarfs. They showed, by means of hydrody- 
namical simulations, that stars belonging to the very high 
helium populations quoted above, found only in the most 
massive clusters, are possibly formed directly from the winds 
of the SAGBs, before mixing with pristin e gas f avour s a less 
extreme chemistry. SAGBs, according to Sicss (2007), have 
helium envelope abundance Y=0.36 - 0.38 after the second 
dredge up, and as such are the only viable candidate for the 
bMSs. Unlike helium, the other yields of SAGBs are the re- 
sult of mass loss and hot bottom burning (HBB) in the phase 
of thermal pulses following the formation of the O-Ne core. 
Common interpretation attributes to them the most extreme 
chemical anomalies found in the clusters having the bMS 
(see, e.g. ID'Ercole et"ai1l20ld ). but a full confirmation will 
have to wait for accurate abundance determinations of stars 
belonging to the bMSs. It is well known that the massive 
AGB yields are dramatically dependent of several uncertain 
input physics parameters, s uch as convection model, mix ing 
and nuclear cross sections l|Ventura fc D 'Antona|[2005aH bl), 
and we can expect similar, or even more extreme, uncer- 
tainties in the SAGB mass range. 

Theoretical models of S AGBs have been first pre - 
sented in a s eries of paper by iGarcia-Berro fc Ibenl (Il994l). 



sented m a s eries or p. 
Ilben et al l (|l997l ). |l 
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1999h . until the recent updates by ISiessI (|20od ) and ISies: 



1999) 
20071 ) 



the carbon burning phase in these models is de- 
scribed in terms of a burning flame, that propogates in- 
wards, until carbon burnin g reaches the center . The main 
effects due to ove rshooting (Gil- Pons et al.ll2007h . and ther- 
mohaline mixing (Siess 20091 ) were also in vestigated and dis- 
cussed. The recent work by [Sicss ( 2010]) presents the first 



1 The range of masses involved depends on the assumptions re- 
garding the overshooting from the convective core during the two 
main phases of nuclear burning: overshooting increases the size of 
the convective core, and thus shifts the masses involved to lower 
values. 



database of yields of SAGBs with different mass and metal- 
licity, and include also the Thermal Pulses (TP) phase that 
follows the e xtinction of carbon b urning. Finally, the in- 
vestigation bv lDohertv et all (2010) outlined the robustness 
of the theoretical description of these evolutionary phases, 
showing that when a sufficiently accurate temporal resolu- 
tion is adopted, results obtained with different codes show 
a sa tisfact o ry agr eement. 

ISiessI (|2010l ) SAGB computations are based on 
physical assumptions ver y different from our own 
l|Ventura fc D'Antonal [2009), that we have already ap- 
plied to the interpretation of GC abundance patterns. Here 
we extend our models to the range of SAGB masses, to 
show how the results depend on detailed physical inputs 
and, mostly, on the assumed mass loss law. In fact, on 
the chemical side, we expect uncertainties on the yields' 
computation, because the mass loss description, the effi- 
ciency of the convective modelling, the choice of the nuclear 
cross-sections, are all expected to play a role in the physical, 
and thus chemical description of the evolutionary phases 
following the formation of the O-Ne core. The fate of these 
stars is also highly uncertain, as it depends on the velocity 
with which the convective mantle is consumed, and, in more 
details, whether this process is c ompleted before electro n 
captures begin inside the core (e.g. lPoelarendes et al.| [2008). 

Our aim is thus to explore the uncertainties associated 
to the yields of these stars, and particularly to the oxy- 
gen and sodium mass fractions, because mostly these two 
elements are investigated in the spectroscopic surveys of 
GC stars. We will also consider the helium content in the 
ejecta, and the mag nesium and aluminum yields, following 
the scheme given i n Ventura fc D 'Antonal |2009) ( but 

sec 

also ID'Ercole et "ail l20ld ). 



2 HOT BOTTOM BURNING: THE 

DIFFERENT FATE OF OXYGEN, SODIUM 
MAGNESIUM AND ALUMINIUM 

Massive AGBs do not obey to the class i c lum inosity vs. core 
mass relationship found by iPaczvnskil (|l970l ). because this 
latter is based on the assumption that a radiative buffer is 
present between the H-burning shell and the bottom of the 
convective envelope, whereas in massive AGBs the bottom 
of the convective envelope eventually overlaps with the H- 
burning shell, so that part of the n uclear energy is generated 
direc tly into the external mantle ^Blocker fc Schoenbernerl 
Il99ll ). This phenomenon, normally indicated as Hot Bot- 
tom Burning, is of paramount importance for the topic of 
the self-enrichment by massive AGBs, because it is an effi- 
cient way of polluting the stellar environment with matter 
that wa s nuclearly processed in repeated proton-capture re- 
actions |Cottrell fc Da Costalll98ll; IVentura et al.ll200ll). 

In a series of papers, IVentura fc D'Anton a ( 2005a. 200^ . 

120091) showed that HBB is expected to operate efficiently in 
all popll AGB models with initial masses M ^ 4Mq, pro- 
vided that a high efficiency convection model is adopted, 
e.g. the Full Spec trum of Turbulence descri ption for tur- 
bulent convection |Canuto fc M azzitclli 199l|). The surface 
chemistry is consequently modified. 

Oxygen is burnt by proton capture; the rate of destruc- 
tion increases for higher temperatures, and is thus stronger 
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in the more massive models, whose core mass is larger. In 
the most massive AGBs, i.e. the high est masses not under- 
going carbon ignition (~ 6 — 6.3M© in lVentura fc D'Antonal 
2008), the final oxygen abundance is a factor of ~ 20 
smaller than the initial m ass fraction (see Figure 2 in 
I Ventura fc D'Antonal IbuOSh . and the average oxygen in the 
ejecta is decreased by a factor ~ 7. When M increases, the 
trend with mass t ends to flatten, and eventua lly is reversed 
(see Table 2 in IVentura fc D'Antonal Iboogl ). because the 
higher temperatures, that tend to decrease the oxygen con- 
tent of the ejecta, are partly compensated by the larger mass 
loss, suffered by the most massive models in the early AGB 
phases, when the surface oxygen was still rather large. This 
is an e ffect of the stee p increase of mass loss with luminosity 
of the lBlockerl l| 19951 ) recipe. 

The behaviour of sodium is more complex, as it is 
a result of a production channel and two destruction re- 
actions, with very different se nsitivity to the temper ature 
in the range of interest here (|Hale et alJ 12001 120041 ). For 
T < 60MK the reaction of proton capture by 22 Ne nu- 
clei dominates, thus producing great amounts of sodium; 
at higher temperatures, the destruction channels via pro- 
ton capture take over, and sodium is destroyed, reaching an 
equilibrium value, at which c reation and burning r a tes ca n- 
cel each other (see Fig. 4 in I Ventura fc D'Antonal l|2008l )). 
Two consequences arise from the above description: 



• Unlike oxygen, the predictions concerning the sodium 
yields are not robust, given the uncertainties of the p-capturc 
reactions of the Ne-Na cycle. 

• The temperature sensitivity of the 23 Na+p reaction, 
much steeper compared to the production channel, if 
confirmed, is necessarily associated to a direct correla- 
tion between the average oxygen and s o dium content 
of the ejecta (iDenissenkoy fc Herwid [20031; iHerwie |2004| ; 



iKarakas fc Lattanziol 120071 ; IVentura fc D'Antonal 12008^ : 
only in the very early AGB phases, when oxygen burning 
is active and the sodium production prevails over destruc- 
tion, oxygen and sodium show an opposite behaviour. 



HBB modifies the surface abundances of magnesium 
and aluminium too. 24 Mg is destroyed when Tb ce ~ 60 — 
70MK, and 27 A1 is produced by a proton-capture chain, that 
also increases the surface mass fracti ons of the two heavy 
magne sium isotopes, 25 Mg and 26 Mg. IVentura fc D'Antonal 
(2008) find that the extent of the Al production (hence, of 
the depletion of magnesium) increases with mass in mas- 
sive AGBs, though an upper limit of [Al/Fe] ~ 1 is reached, 
which is again due to the high mass loss experienced by 
the most massive stars at the beginning of the AGB phase. 
These results were obtained by adopting the upper limits for 
the two proton capture reactions by the he avy magnesium 
isotop es, although the recent analysis by iD'Ercole et al.l 
(|2010r i on M4 recommends the use of the standard cross 
sections. The uncertainties associated to the relevant cross- 
sections of t he Mg-Al chain rende r the r esults uncertain by 
~ 0.3 dex (IVentura fc D'Antonal (120081 1. but see also the 
discussion in Izzard et al. 1 20071 )'). 



3 THE MODELS 

3.1 Physical and chemical inputs 

To investigate the nucleosynthesis at the bottom of the ex- 
ternal mantle of SAGB stars after the carbon burning phase, 
we calculate models with the standard chemistry typical of 
intermediate metallicity GCs, i.e. Z=0.0 01, Y=0.24, and an 
q-enh anced mixture, with [a/Fe]=+0.4 l|Grevesse fc Sauvall 
Il998h . 

The nuclear network includes the most important re- 
actions involving all the chemical species up to silicon. 
Most of the cross-sections were taken from the NACRE 
com pilation (jAngulo et all 19991); t o mai ntain consistency 
with IVentura fc D'Antonal (|2005al lbl. |200S| ) we used the up- 
per limits for the proton capture reactions by 25 Mg and 
26 Mg. This choice maximizes the aluminum production; a 
full description of the Mg-Al production and destruction 
will be given in a forthcoming paper. The reaction rates of 
the 14 N(p,7) 15 Q reaction were take n from iFormicola et alJ 
J2004h. whereas the cross sections bv lHale et all (j2002f ) and 
lHale et al.l (|2004f ) were adopted for the reactions of the Ne- 
Na cycle. Sodium production was maximized by adopting 
the upper limit for the proton capture reaction by 22 Ne nu- 
clei. 

Mass loss was modelled according to iBlockerl |l995), 
while some compa ri son s equences follow the recipe by 
IVassiliadis fc Wood! (|l993l ). Nuclear burning and convec- 
tive mixing were couple d by means of a diffusive approach 
dCloutman fc Eolllli976). Convective eddies are allowed to 
overshoot into regions of radiative stability during both hy- 
drogen and helium nuclear burning in the core; this extra- 
mixing is simulated by an exponential decay of velocities be- 
yond the formal border fixed via the classic Schwartzschild 
criterium, with an e-folding distance C,H P , with ( — 0.02, 
in agreement with a calibration based on the observed 
w idths of the main se quences of open clusters, presented 
in lVentura et al.l (| 19981 ). The range of masses investigated is 
6.5-8 solar masses; this range is to be compare d with the in - 
terval 8-9M© found for the same chemistry bv lSiess! (|2010h . 
The different assumptions on overshooting are the cause of 
the difference in the mass range of SAGBs. The involved 
core masses in fact are very similar in the two investigations 
(see Figure [5]). The results are also in agre ement with the 
explor ation on the effects of overshooting by I Gil- Pons et all 

<Eqq3). 



3.2 The Carbon burning phase 

The main physical properties concerning the C-burning 
phase are presented in Tab. [TJ where for the various masses 
involved we show the time elapsed from the end of the core- 
He burning phase to the ignition of the first C-burning 
episode, the mass of the CO core when carbon burning 
starts, the mass of the layer at which carbon is ignited, 
the maximum luminosity produced, and the innermost layer 
reached by the bottom of the surface convective zone during 
the second dredge-up. 

All the models undergo an off-center carbon ignition; 
the point at which the nuclear energy release is highest is 
more internal the higher is the mass, ranging from ~ O.6A/0 
for M=6.5M Q , down to ~ 0.2M© for M=8M . This early 
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Table 1. Carbon burning properties 



M/Mq 


TCB-HeB (10 5 )yr 


M co 




Lc 


M SDU 


6.5 


2.3 


1.06 


0.60 


7.1(6) 


1.15 


7.0 


2.0 


1.12 


0.50 


2.3(6) 


1.20 


7.5 


1.7 


1.18 


0.36 


1.0(6) 


1.30 


8.0 


1.4 


1.29 


0.21 


6.3(5) 


1.40 




0.4 



1.85x10" 1.9x10° 
Time 



Figure 1. Kippenhan diagram showing the temporal behaviour 
of the development of convective regions during the C-burning 
phase of a 7Mq model. 



phase of carbon burning is followed by the development of 
a convective flame, that proceeds inwards. In all the mod- 
els investigated, this instability region eventually reaches the 
center, and forms a O-Ne core, with the only exception of the 
6.5Mq model, that is left with a core made up of carbon and 
oxygen; in this latter model carbon burning is aborted, in 
analogy to the behaviour of the masses just above the limit 
for carbon ignition, described bv lDohertv et all (2010). The 
formation of the O-Ne core is followed by other off-center 
carbon burning episodes, whose intensity becomes progres- 
sively lower, after which the TPs phase begins. Fig. [1] shows 
a Kippenhahn diagram for the model with mass 7Mq, in- 
dicating the borders of the convective shells that form as a 
consequence of carbon ignition, and the base of the exter- 
nal convective mantle. A detailed analysis of the C-burning 
phase is beyond the scop es of the pre sent work: the inter- 
ested reader may find in ISiessI (2006) a detailed physical 
description of the C-burning phase in super-AGB stars. 



3.3 How much helium is produced by SAGBs? 

The second dredge-up is a common feature of all the models 
investigated; the evolutionary stage at which it takes place 
depends on the initial mass: it is prior to C-burning in all the 
models, apart form the 8M© case, in which it occurs after the 
formation of the ONe core. This behaviour is in agreement 
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Figure 2. The heli um content of the ejecta of the massive AGB 
(open squares, from IVentura fc D'Antonal20 09') and of the SAGB 
models calculated in this work (full squares) is shown as a function 
of the core mass. The numbers next to the points indicat e the 
initial mass of the model. The values for the SAGBs by ISiessI 
(2010) (full triangles) for similar chemistry are shown. The full 
dot indicates how the helium content from the 7M0 decreases 
when a shallower mass loss formulation is adopted. 



with the investigation bv lGarcia- Berro et all i| 19971 ). On the 
chemical side, the main consequence of the second dredge-up 
is the increase in the helium mass fraction, up to Y — 0.36 
(see column 4 in Table [2J). This is the maximum enrich- 
ment of helium that is produced by SAGBs according to the 
present investigation. It is in reson able a g reem ent, though a 
bit lower, than the helium yields bv lSiess! ( 2010l ) for the same 
metallicity. Figure [2] sho ws the abundances in the ej ecta of 
the massive AGBs from IVentura fc D'Antonal ((2009) (open 
squares) and of the present models (full squares) , as a func- 
tion of the core mass (CO cores from 4 to 6.5M0 , ONe cores 
for larger mass). 

Our SAGB models do not provide helium abundances 
as larg e as Y^0.40 invok ed to interpret the bMS in 
to Ce n (|Sollima et all 120071 ') and NGC 2808 |Piotto et all 
l2007h . Notice however that the values depend on the inter- 
pretat ion of the MS colors , on which uncertainties may be 
large |Portinari et al.ll2010l i. 
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Table 2. Properties of SAGB models and average abundances in the ejecta 



M/M 


M c /M 


T max / 1Q 8 
bcc / 


Y 


[ ie O/Fe] 


[Na/Fe] 


[Mg/Fe] 


[ 27 A1/Fe] 


A(Li) 


M law 


Mf ot /M 


m£/m 


NTP 


6.5 


1.08 


1.16 


0.352 


-0.24 


0.32 


0.226 


0.80 


2.36 


B95 


1.70 


1.10 


32 


7.0 


1.20 


1.20 


0.358 


-0.15 


0.39 


0.253 


0.74 


2.12 


B95 


1.99 


1.23 


38 


7.5 


1.27 


1.27 


0.359 


0.01 


0.67 


0.289 


0.57 


2.75 


B95 


2.15 


1.29 


48 


8.0 


1.36 


1.47 


0.344 


0.20 


1.00 


0.290 


0.40 


4.39 


B95 


2.20 


1.38 


53 


7.0 


1.20 


1.34 


0.350 


-0.69 


0.17 


0.09 


0.54 


2.58 


VW93 


2.15 


1.24 


348 




Figure 3. Left panel: The variation during the evolution of the sodium abundance at the surface of models with initial masses 6M 
(light solid line), 6.5M (dotted), 7M (dashed), 7.5M (dot-dashed), 8M (solid); Right panel: the same as the left panel, but for the 
oxygen surface mass fraction. 



3.4 The effects of HBB 

Tab. [2] shows the main properties of the TP phase of our 
models, i.e. the core mass at the first thermal pulse, the 
maximum temperature reached at the base of the convective 
mantle, the average chemistry of the ejecta, the mass of the 
core and of the whole star when the computations stopped, 
and the number of thermal pulses experienced. 

During the TP phase no extra-mixing was consid- 
ered from any convective border, and no appreciable third 
dredge- up was experienced by any of our models: the surface 
chemistry changes only for the effects of HBB. 

Fig. [3] shows the behaviour of the surface oxygen 
(right) and sodium (left), as the envelope mass is con- 
sumed. To show the continuity with the previous investiga- 
tions, limited to models not undergoing any carbon burn- 
ing, we also show the evolution of a 6M model, taken 
from lVentura fc D'Antonal (|2008l ). Oxygen is destroyed in all 
models during the thermal pulses phase. Despite the higher 
temperatures attained at the bottom of the convective en- 
velope in the more massive models (see Table (2}, the higher 
mass loss experienced (due to the larger luminosity) favours 
a flatter oxygen vs. mass relationship, because a lot of mass 
is lost before the oxygen is burnt in great quantities. 

The average oxygen content of the ejecta, which in Ta- 



ble [2] is expressed in terms of [0/Fej3 is found to increase 
with mass; this trend agree s with the previous findings by 
I Ventura fc D'Antonal (|2009l ) (most massive AGBs) and by 
bicss (SAGBs), and indicates that the stars whose 

ejecta show the most extreme chemistries are those close to 
the limit for carbon ignition. 

The behaviour of sodium is more tricky. The maxi- 
mum sodium produced increases with mass (as a conse- 
quence of the larger temperatures attained at the bottom 
of the convective envelope), with the only exception of the 
8M model, that achieves the second dredge-up after the 
C-burning phase, and evolves rapidly to high temperatures: 
in this case sodium begins to be destroyed before the early 
increase, normally associated to the conversion of the 22 Ne 
transported to the surface by the second dredge-up. Simi- 
larly to oxygen, the sodium content of the ejecta increases 
with mass (see Table [2]). 

The simultaneous behaviour of oxygen and sodium 
abundance at the surface of these stars along their evolu- 
tion is shown in Fig(4] We used the quantity [X/Fe] on both 
axis. We see that all the models follow the same qualitative 

2 Here we use the standard notation, by indicating with [X/Fe] 
the quantity log(X/Fe) - \og(X/Fe) Q . 
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Figure 4. The simultaneous variation of the surface abundances 
of oxygen and sodium in the SAGB models presented in this in- 
vestigation. The labels are the same as in the two panels of Fig[3] 
For clarity reasons, the path followed by the 8Mq model is indi- 
cated by solid triangles, and is terminated by the full square 



Figure 5. The simultaneous variation of the surface abundances 
of magnesium and aluminium in the SAGB models presented in 
this investigation. The labels are the same as in the two panels 
of Fig|3] 



behaviour, with an initial anticorrelated phase, during which 
oxygen is reduced whereas sodium increases, and a following 
phase, when both elements are destroyed at the bottom of 
the external envelope. 

Oxygen is destroyed with continuity until the envelope 
is removed completely, whereas sodium follows an approxi- 
mately asymptotic behaviour, that corresponds to the afore 
mentioned balance between the production and the destruc- 
tion channels. The curves corresponding to the more massive 
models (see in particular the evolution of the 8Mq model, 
indicated by solid triangles and by the full square pointing 
the end of the evolution), stop at a stage when the oxygen 
and sodium abundances are still large, because the strong 
mass loss favoured a rapid consumption of the whole en- 
velope, before the surface mass fractions of both elements 
could diminish to nominal values. 

The variation of the surface contents of magnesium and 
27 A1 is shown in FigO We note that a stronger production 
of aluminium is achieved in the models of smaller mass, be- 
cause in the most massive stars the strong mass loss prevents 
an advanced Mg-Al nucleosynthesis, and only a modest in- 
crease in the surface aluminium is reached. 



4 HOW SENSITIVE ARE THE YIELDS TO 
THE DETAILS OF SAGB MODELLING? 

Before turning to the impact of these results on the interpre- 
tation of the chemical anomalies observed in GC stars, and 
on the self-enrichment scenario by massive AGBs, we try 
to understand how robust these results are, and their sen- 
sitivity to the various uncertainties (convection, mass loss, 

cross-sections) affecting the SAG B modelling. 

The recent exploration by iDohertv et alj (|2010l ) re- 
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Figure 6. The average [ O/Fe] and [N a/Fel content of the ejecta 
of the SAGB models by ISiessi ll2010h (full triangles), and those 
calculated in this work (ful l squares ). Open squares denote the 
massive AGB yields bv lVentura &: D'Antonal &200$) . The num- 
bers next to the points indicate the initial mass of the model, 
and the core mass (in parenthesis). The arrow indicates how the 
yields of the 7Mq model change when the lower mass loss formu- 
lation is adopted (see Table 1). 
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vealed a great homogeneity among the results obtained with 
different codes, that confirm the results concerning the C- 
burning phase already present in the literature. These mod- 
els are limited to the evolutionary phases preceeding the first 
thermal pulse, whereas as far as we know the only full inves- 
tigation extend ed to the whole SAGB evolution was given 
bv lSiess! (|2010h . 

Fig. [6] shows a co mparison between the yields of SAGB 
models by Siesd (j2010l ) (indicated as full triangles) and those 
found in the present investigation (open squares), in the 
[Na/Fe] vs [O/Fe] plane. In both cases the points move 
towards the to p-righ t of th e plane with increasing mass. 
The yields by Siesj l)2010l ) show lower contents of oxy- 
gen and sodium. The differences in [O/Fe] are mostly due 
to the different mo d elling of mass lo ss. The trea tment by 
IVassiliadis fc Wood! |l993h . used by ISiessl (|20ld ). favours 
much smaller rates, thus leading to a very advanced nu- 
cleosynthesis at the bott om of the oute r convective zone. 
Conversely, the recipe by iBlockerl (|l995l ) determines much 
higher Ms, so that the mass of the envelope is lost before 
a great destruction of oxygen can be achieved. This is con- 
firmed by the full dot in Fig El that indicates the yield found 
in the TM© model when the Vassiliadis fc Wood! (|l993h rate 
is adopted: the theoretical poi nt is shifted in the range of 
oxygen abundances predicted by Sicss (2010). The remaining 
difference is due to the fact that our initial mixture is alpha- 
enhanced w ith [a/Fe] = +0.4, whereas the mixture used by 
ISiessl (|2010T ) is solar-scaled. The sodium yields remain dif- 
ferent, our [Na/Fe] being sistematically larger. This can be 
attributed to the different nuclear reaction rates adopted, in 
particular for what concerns the proton capture reaction by 
22 Ne nucle i, for w hich we adopted the highest rates from the 
lHale et al.l (|2002T) compilatio ns, whereas the recommended 



rate was used bv lSiessI 

The comparison for what concerns the magnesium and 
aluminium content is more tricky, because of the differences 
in the assumptions made, that involve the nuclear cross- 
sections, the initial magnesium abundances, and the mass 
loss description. The Al yields bv lSiessI (|2010f ) decrease with 
mass, starting form [Al/Fe]=0.27 for the M=8M model, 
down to [Al/Fe]=0.18 for M=9M Q . Our yields show the 
same trend (see col. 8 in Table [2]), but the values are sistem- 
atically higher, mainly because our alpha-enhanced mixture 
has a higher initial content of magnesium. 

This investigation indicates that the oxygen and sodium 
yields in SAGBs vary considerably according to the treat- 
ment of mass loss: when a modest mass loss is adopted, 
extreme chemistries are possible, because a very advanced 
nucleosynthesis can be achieved at the bottom of the con- 
vective envelope. Convective modelling is also expected to 
play a role in de t ermin ing the yields of SAGBs (see the dis- 
cussion m ISiessl l|2010l) ). although the effects are less im- 
portant than for AGBs, because a more efficient convective 
modelling would determine higher temperatures, but also a 
larger mass loss, and the two effects tend to compensate for 
what concerns the stage of nucleosynthesis achieved. 

In the near future, it could be possible to calibrate the 
mass loss rate by achieving constraints from the abundances 
of stars belonging to the bMS in the clusters uj Cen and 
NGC 2808, in the hypothesis that the bMS stars are di- 
rectly formed from t he ashes of SAGBs l)Pumo et al.ll20oljl : 
iD'Ercole etailbOOSl ). 



5 CONCLUSIONS 

We present the evolution of the main physical and chemical 
properties of stars with initial mass in the range 6.5-8Mq, 
which, according to our choice for the overshooting from 
the convective core during the H- and He-burning phases, 
undergo carbon ignition in conditions of partial degener- 
acy. Our findings concerning the modalities of carbon burn- 
ing, the inwards propagation of the convective flame that 
is developed, and the final chemistry of the O-Ne core, are 
in agreement with the results presented by other research 
groups. 

The same homogeneity cannot be expected for the fol- 
lowing TP phase, because its physical and chemical evo- 
lution is determined by the treatment of two still poorly 
known physical processes: mass loss and turbulent convec- 
tion. We focused our analysis on the chemistry of the ejecta 
of SAGBs, and in particular on the oxygen, sodium and alu- 
minum content, that have been shown to vary among stars 
in globular clusters. 

The models presented here attain extremely large tem- 
peratures at the bottom of the convective envelope, which, 
in turn, favour an advanced nucleosynthesis, i.e. very small 
abundances of oxygen and sodium, and a strong production 
of aluminum. These findings, whic h are i n fairl y good agree- 
ment with the results obtained bv lSiessI (|2010T). hold in case 
that m ass loss is modelled ac cording t o Vassiliadis fc Wood! 
(|l993h . Conversely, when the lBlockerl (|l995l ) recipe is used, 
the nucleosynthesis at the bottom of the external mantle is 
halted by the general cooling of the structure favoured by 
the fast consumption of the envelope, so that a less extreme 
chemistry is achieved. In this latter case higher abundances 
of oxygen and sodium are produced, and, in the overall con- 
text of a possible pollution in GCs produced by massive 
AGBs and super- AGBs, gas with the most extreme chem- 
istry is expected to be ejected by masses around 6Mq . 

Helium is produced in great quantities in all cases, with 
a maximum of Y~ 0.36 for the highest masses. In the frame- 
work of our models, abundances larger than this value for 
the extreme populations in the most massive clusters cannot 
be explained by formation from SAGB ashes. However, the 
values derived from o bservations are still largely uncertain 
jPortinari et al.ll2010l ). 

If the bMS stars in massive GCs are formed directly 
from the ejecta of SAGBs, future spectroscopic analysis will 
help understanding the reliability of the self-enrichment sce- 
nario by massive AGBs, and, eventually, will help to cali- 
brate the poorly constrained mass loss in the SAGB phase. 
In the end, this calibration may also be important to con- 
strain the mass range in which we expect that the core grows 
up to begin electron capture on the Ne nuclei. Indirect infor- 
mations on the formation of neutron stars by the e-capture 
supernova channel in sin gle stars will be eventually achieved 
jPoelarendes et al.ll200sh . 
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